The gut microbiome reflects ancestry despite dietary shifts across a hybrid zone.

The microbiome is critical to an organism's phenotype, and its composition is shaped by, and a driver of, eco-evolutionary interactions. We investigated how host ancestry, habitat and diet shape gut microbial composition in a mammalian hybrid zone between Neotoma lepida and N. bryanti that occurs across an ecotone between distinct vegetation communities. We found that habitat is the primary determinant of diet, while host genotype is the primary determinant of the gut microbiome-a finding further supported by intermediate microbiome composition in first-generation hybrids. Despite these distinct primary drivers, microbial richness was correlated with diet richness, and individuals that maintained higher dietary richness had greater gut microbial community stability. Both relationships were stronger in the relative dietary generalist of the two parental species. Our findings show that host ancestry interacts with dietary habits to shape the microbiome, ultimately resulting in the phenotypic plasticity that host-microbial interactions allow.

Phylogenomic analysis of Syngnathidae reveals novel relationships, origins of endemic diversity and variable diversification rates.

BACKGROUND: Seahorses, seadragons, pygmy pipehorses, and pipefishes (Syngnathidae, Syngnathiformes) are among the most recognizable groups of fishes because of their derived morphology, unusual life history, and worldwide distribution. Despite previous phylogenetic studies and recent new species descriptions of syngnathids, the evolutionary relationships among several major groups within this family remain unresolved. RESULTS: Here, we provide a reconstruction of syngnathid phylogeny based on genome-wide sampling of 1314 ultraconserved elements (UCEs) and expanded taxon sampling to assess the current taxonomy and as a basis for macroevolutionary insights. We sequenced a total of 244 new specimens across 117 species and combined with published UCE data for a total of 183 species of Syngnathidae, about 62% of the described species diversity, to compile the most data-rich phylogeny to date. We estimated divergence times using 14 syngnathiform fossils, including nine fossils with newly proposed phylogenetic affinities, to better characterize current and historical biogeographical patterns, and to reconstruct diversification through time. We present a phylogenetic hypothesis that is well-supported and provides several notable insights into syngnathid evolution. We found nine non-monophyletic genera, evidence for seven cryptic species, five potentially invalid synonyms, and identified a novel sister group to the seahorses, the Indo-Pacific pipefishes Halicampus macrorhynchus and H. punctatus. In addition, the morphologically distinct southwest Pacific seahorse Hippocampus jugumus was recovered as the sister to all other non-pygmy seahorses. As found in many other groups, a high proportion of syngnathid lineages appear to have originated in the Central Indo-Pacific and subsequently dispersed to adjoining regions. Conversely, we also found an unusually high subsequent return of lineages from southern Australasia to the Central Indo-Pacific. Diversification rates rose abruptly during the Middle Miocene Climate Transition and peaked after the closure of the Tethys Sea. CONCLUSIONS: Our results reveal a previously underappreciated diversity of syngnathid lineages. The observed biogeographic patterns suggest a significant role of the southern Australasian region as a source and sink of lineages. Shifts in diversification rates imply possible links to declining global temperatures, the separation of the Atlantic and Pacific faunas, and the environmental changes associated with these events.

Phylogeny of the supertribe Nebriitae (Coleoptera, Carabidae) based on analyses of DNA sequence data.

The phylogeny of the carabid beetle supertribe Nebriitae is inferred from analyses of DNA sequence data from eight gene fragments including one nuclear ribosomal gene (28S), four nuclear-protein coding genes (CAD, topoisomerase 1, PEPCK, and wingless), and three mitochondrial gene fragments (16S + tRNA-Leu + ND1, COI ("barcode" region) and COI ("Pat/Jer" region)). Our taxon sample included 264 exemplars representing 241 species and subspecies (25% of the known nebriite fauna), 39 of 41 currently accepted genera and subgenera (all except Notiokasis and Archileistobrius), and eight outgroup taxa. Separate maximum likelihood (ML) analyses of individual genes, combined ML analyses of nuclear, nuclear protein-coding, and mitochondrial genes, and combined ML and Bayesian analyses of the eight-gene-fragment matrix resulted in a well-resolved phylogeny of the supertribe, with most nodes in the tree strongly supported. Within Nebriitae, 167 internal nodes of the tree (out of the maximum possible 255) are supported by maximum-likelihood bootstrap values of 90% or more. The tribes Notiophilini, Opisthiini, Pelophilini, and Nebriini are well supported as monophyletic but relationships among these are not well resolved. Nippononebria is a distinct genus more closely related to Leistus than Nebria. Archastes, Oreonebria, Spelaeonebria, and Eurynebria, previously treated as distinct genera by some authors, are all nested within a monophyletic genus Nebria. Within Nebria, four major clades are recognized: (1) the Oreonebria Series, including eight subgenera arrayed in two subgeneric complexes (the Eonebria and Oreonebria Complexes); (2) the Nebriola Series, including only subgenus Nebriola; (3) the Nebria Series, including ten subgenera arrayed in two subgeneric complexes, the Boreonebria and Nebria Complexes, with the latter further subdivided into three subgeneric subcomplexes (the Nebria, Epinebriola, and Eunebria Subcomplexes)); and (4) the Catonebria Series, including seven subgenera arrayed in two subgeneric complexes (the Reductonebria and Catonebria Complexes). A strong concordance of biogeography with the inferred phylogeny is noted and some evident vicariance patterns are highlighted. A revised classification, mainly within the Nebriini, is proposed to reflect the inferred phylogeny. Three genus-group taxa (Nippononebria, Vancouveria and Archastes) are given revised status and seven are recognized as new synonymies (Nebriorites Jeannel, 1941 and Marggia Huber, 2014 = Oreonebria Daniel, 1903; Pseudonebriola Ledoux & Roux, 1989 = Boreonebria Jeannel, 1937; Patrobonebria Bänninger, 1923, Paranebria Jeannel, 1937 and Barbonebriola Huber & Schmidt, 2017 = Epinebriola Daniel & Daniel, 1904; and Asionebria Shilenkov, 1982 = Psilonebria Andrewes, 1923). Six new subgenera are proposed and described for newly recognized clades: Parepinebriola Kavanaugh subgen. nov. (type species: Nebria delicata Huber & Schmidt, 2017), Insulanebria Kavanaugh subgen. nov. (type species: Nebria carbonaria Eschscholtz, 1829), Erwinebria Kavanaugh subgen. nov. (type species Nebria sahlbergii Fischer von Waldheim, 1828), Nivalonebria Kavanaugh subgen. nov. (type species: Nebria paradisi Darlington, 1931), Neaptenonebria Kavanaugh subgen. nov. (type species: Nebria ovipennis LeConte, 1878), and Palaptenonebria Kavanaugh subgen. nov. (type species: Nebria mellyi Gebler, 1847). Future efforts to better understand relationships within the supertribe should aim to expand the taxon sampling of DNA sequence data, particularly within subgenera Leistus and Evanoleistus of genus Leistus and the Nebria Complex of genus Nebria.

Genome-wide analyses reveal drivers of penguin diversification.

Penguins are the only extant family of flightless diving birds. They currently comprise at least 18 species, distributed from polar to tropical environments in the Southern Hemisphere. The history of their diversification and adaptation to these diverse environments remains controversial. We used 22 new genomes from 18 penguin species to reconstruct the order, timing, and location of their diversification, to track changes in their thermal niches through time, and to test for associated adaptation across the genome. Our results indicate that the penguin crown-group originated during the Miocene in New Zealand and Australia, not in Antarctica as previously thought, and that Aptenodytes is the sister group to all other extant penguin species. We show that lineage diversification in penguins was largely driven by changing climatic conditions and by the opening of the Drake Passage and associated intensification of the Antarctic Circumpolar Current (ACC). Penguin species have introgressed throughout much of their evolutionary history, following the direction of the ACC, which might have promoted dispersal and admixture. Changes in thermal niches were accompanied by adaptations in genes that govern thermoregulation and oxygen metabolism. Estimates of ancestral effective population sizes (Ne ) confirm that penguins are sensitive to climate shifts, as represented by three different demographic trajectories in deeper time, the most common (in 11 of 18 penguin species) being an increased Ne between 40 and 70 kya, followed by a precipitous decline during the Last Glacial Maximum. The latter effect is most likely a consequence of the overall decline in marine productivity following the last glaciation.

Distinct patterns of hybridization across a suture zone in a coral reef fish (Dascyllus trimaculatus).

Hybrid zones are natural laboratories for investigating the dynamics of gene flow, reproductive isolation, and speciation. A predominant marine hybrid (or suture) zone encompasses Christmas Island (CHR) and Cocos (Keeling) Islands (CKE), where 15 different instances of interbreeding between closely related species from Indian and Pacific Oceans have been documented. Here, we report a case of hybridization between genetically differentiated Pacific and Indian Ocean lineages of the three-spot dascyllus, Dascyllus trimaculatus (Rüppell, 1829). Field observations indicate there are subtle color differences between Pacific and Indian Ocean lineages. Most importantly, population densities of color morphs and genetic analyses (mitochondrial DNA and SNPs obtained via RADSeq) suggest that the pattern of hybridization within the suture zone is not homogeneous. At CHR, both color morphs were present, mitochondrial haplotypes of both lineages were observed, and SNP analyses revealed both pure and hybrid genotypes. Meanwhile, in CKE, the Indian Ocean color morphs were prevalent, only Indian Ocean mitochondrial haplotypes were observed, and SNP analysis showed hybrid individuals with a large proportion (~80%) of their genotypes assigning to the Indian Ocean lineage. We conclude that CHR populations are currently receiving an influx of individuals from both ocean basins, with a greater influence from the Pacific Ocean. In contrast, geographically isolated CKE populations appear to be self-recruiting and with more influx of individuals from the Indian Ocean. Our research highlights how patterns of hybridization can be different at scales of hundreds of kilometers, due to geographic isolation and the history of interbreeding between lineages.

Patterns of Genomic Divergence and Signals of Selection in Sympatric and Allopatric Northeastern Pacific and Sea of Cortez Populations of the Sargo (Anisotremus davidsonii) and Longjaw Mudsucker (Gillichthys mirabilis).

Studying how isolation can impact population divergence and adaptation in co-distributed species can bring us closer to understanding how landscapes affect biodiversity. The Sargo, Anisotremus davidsonii (Haemulidae), and the Longjaw mudsucker, Gillichthys mirabilis (Gobiidae), offer a notable framework to study such mechanisms as their Pacific populations cross phylogeographic breaks at Point Conception, California, United States, and Punta Eugenia, Mexico, and are separated to those in the Sea of Cortez by the Baja California peninsula. Here, thousands of loci are genotyped from 48 Sargos and 73 mudsuckers using RADseq to characterize overall genomic divergence, and search for common patterns of putatively neutral and non-neutral structure based on outlier loci among populations with hypothesized different levels of isolation. We further search for parallels between population divergence and the total proportion of outliers, outlier FST distribution, and the proportion of outliers matching coding regions in GenBank. Statistically significant differentiation is seen across Point Conception in mudsucker (FST = 0.15), Punta Eugenia in Sargo (FST = 0.02), and on either side of the Baja California peninsula in both species (FST = 0.11 and 0.23, in Sargo and mudsucker, respectively). Each species shows structure using neutral and non-neutral loci. Finally, higher population divergence yields a more even distribution of outliers along their differentiation range but does not always translate into higher outlier proportions or higher rates in which outliers are matched to coding regions. If repeated in similar systems, observed genomic patterns might reveal speciation signatures in diverse networks of population isolation.

Genetic variation and admixture of red-eared sliders (Trachemys scripta elegans) in the USA.

The most ubiquitous, abundant, and invasive turtle on Earth, Trachemys scripta elegans (TSE, "red-eared slider"), is one of four taxa in a clade that is native to the USA and adjacent Mexico (three subspecies of Trachemys scripta plus Trachemys gaigeae). The present range-wide study of this clade is based on 173 known-locality mtDNA sequences combined with ddRAD libraries for 43 samples emphasizing the western part of the range of TSE, its contact with that of T. gaigeae, and anthropogenic hybrids between TSE and T. s. scripta. The data presented here are the first to sample the TSE × T. s. scripta intergrade zone or TSE × T. s. scripta crosses from introduced turtles. In the western part of its range (New Mexico and Texas), most samples of TSE from the Pecos River have mtDNA haplotypes matching T. gaigeae. Structure analysis of SNPs from the ddRAD show evidence of genetic admixture between T. gaigeae and TSE in all included samples from the Rio Grande and Pecos River. These populations also exhibit T. gaigeae-like head stripes, i.e., a postorbital marking that does not reach the eye. The genetic and morphological data are thereby reconciled, as both suggest that these TSE are intergrades. We recommend that these populations continue to be considered TSE, despite the admixture with T. gaigeae. In the Eastern United States, some samples of the morphologically intermediate subspecies T. s. troostii are not genetically distinct from TSE and some samples share morphological characters and genetic affinities with T. s. scripta. Based on these observations we conclude that the taxon T. s. troostii represents intergrades between TSE and T. s. scripta and should not be considered a valid taxon. Near the already established part of the intergrade zone between TSE and T. s. scripta, TSE mtDNA haplotypes have naturally introgressed into typical-looking samples of T. s. scripta in Georgia. Hybrids between introduced TSE and T. s. scripta are also confirmed deeper within the natural range of T. s. scripta in South Carolina and Virginia. Given the examples of feral hybrids deep within its range shown here and elsewhere, the threat of genetic pollution of T. s. scripta by feral TSE is established.

Whole-genome assembly of the coral reef Pearlscale Pygmy Angelfish (Centropyge vrolikii).

The diversity of DNA sequencing methods and algorithms for genome assemblies presents scientists with a bewildering array of choices. Here, we construct and compare eight candidate assemblies combining overlapping shotgun read data, mate-pair and Chicago libraries and four different genome assemblers to produce a high-quality draft genome of the iconic coral reef Pearlscale Pygmy Angelfish, Centropyge vrolikii (family Pomacanthidae). The best candidate assembly combined all four data types and had a scaffold N50 127.5 times higher than the candidate assembly obtained from shotgun data only. Our best candidate assembly had a scaffold N50 of 8.97 Mb, contig N50 of 189,827, and 97.4% complete for BUSCO v2 (Actinopterygii set) and 95.6% complete for CEGMA matches. These contiguity and accuracy scores are higher than those of any other fish assembly released to date that did not apply linkage map information, including those based on more expensive long-read sequencing data. Our analysis of how different data types improve assembly quality will help others choose the most appropriate de novo genome sequencing strategy based on resources and target applications. Furthermore, the draft genome of the Pearlscale Pygmy angelfish will play an important role in future studies of coral reef fish evolution, diversity and conservation.

Ancient DNA Resolves the History of Tetragnatha (Araneae, Tetragnathidae) Spiders on Rapa Nui.

Rapa Nui is one of the most remote islands in the world. As a young island, its biota is a consequence of both natural dispersals over the last ~1 million years and recent human introductions. It therefore provides an opportunity to study a unique community assemblage. Here, we extract DNA from museum-preserved and newly field-collected spiders from the genus Tetragnatha to explore their history on Rapa Nui. Using an optimized protocol to recover ancient DNA from museum-preserved spiders, we sequence and assemble partial mitochondrial genomes from nine Tetragnatha species, two of which were found on Rapa Nui, and estimate the evolutionary relationships between these and other Tetragnatha species. Our phylogeny shows that the two Rapa Nui species are not closely related. One, the possibly extinct, T. paschae, is nested within a circumtropical species complex (T. nitens), and the other (Tetragnatha sp. Rapa Nui) appears to be a recent human introduction. Our results highlight the power of ancient DNA approaches in identifying cryptic and rare species, which can contribute to our understanding of the global distribution of biodiversity in all taxonomic lineages.

Sequence-Based Analysis of Thermal Adaptation and Protein Energy Landscapes in an Invasive Blue Mussel (Mytilus galloprovincialis).

Adaptive responses to thermal stress in poikilotherms plays an important role in determining competitive ability and species distributions. Amino acid substitutions that affect protein stability and modify the thermal optima of orthologous proteins may be particularly important in this context. Here, we examine a set of 2,770 protein-coding genes to determine if proteins in a highly invasive heat tolerant blue mussel (Mytilus galloprovincialis) contain signals of adaptive increases in protein stability relative to orthologs in a more cold tolerant M. trossulus. Such thermal adaptations might help to explain, mechanistically, the success with which the invasive marine mussel M. galloprovincialis has displaced native species in contact zones in the eastern (California) and western (Japan) Pacific. We tested for stabilizing amino acid substitutions in warm tolerant M. galloprovincialis relative to cold tolerant M. trossulus with a generalized linear model that compares in silico estimates of recent changes in protein stability among closely related congeners. Fixed substitutions in M. galloprovincialis were 3,180.0 calories per mol per substitution more stabilizing at genes with both elevated dN/dS ratios and transcriptional responses to heat stress, and 705.8 calories per mol per substitution more stabilizing across all 2,770 loci investigated. Amino acid substitutions concentrated in a small number of genes were more stabilizing in M. galloprovincialis compared with cold tolerant M. trossulus. We also tested for, but did not find, enrichment of a priori GO terms in genes with elevated dN/dS ratios in M. galloprovincialis. This might indicate that selection for thermodynamic stability is generic across all lineages, and suggests that the high change in estimated protein stability that we observed in M. galloprovincialis is driven by selection for extra stabilizing substitutions, rather than by higher incidence of selection in a greater number of genes in this lineage. Nonetheless, our finding of more stabilizing amino acid changes in the warm adapted lineage is important because it suggests that adaption for thermal stability has contributed to M. galloprovincialis' superior tolerance to heat stress, and that pairing tests for positive selection and tests for transcriptional response to heat stress can identify candidates of protein stability adaptation.

Introgression and selection shaped the evolutionary history of sympatric sister-species of coral reef fishes (genus: Haemulon).

Closely related marine species with large overlapping ranges provide opportunities to study mechanisms of speciation, particularly when there is evidence of gene flow between such lineages. Here, we focus on a case of hybridization between the sympatric sister-species Haemulon maculicauda and H. flaviguttatum, using Sanger sequencing of mitochondrial and nuclear loci, as well as 2422 single nucleotide polymorphisms (SNPs) obtained via restriction site-associated DNA sequencing (RADSeq). Mitochondrial markers revealed a shared haplotype for COI and low divergence for CytB and CR between the sister-species. On the other hand, complete lineage sorting was observed at the nuclear loci and most of the SNPs. Under neutral expectations, the smaller effective population size of mtDNA should lead to fixation of mutations faster than nDNA. Thus, these results suggest that hybridization in the recent past (0.174-0.263 Ma) led to introgression of the mtDNA, with little effect on the nuclear genome. Analyses of the SNP data revealed 28 loci potentially under divergent selection between the two species. The combination of mtDNA introgression and limited nuclear DNA introgression provides a mechanism for the evolution of independent lineages despite recurrent hybridization events. This study adds to the growing body of research that exemplifies how genetic divergence can be maintained in the presence of gene flow between closely related species.

Global divergence of the human follicle mite Demodex folliculorum: Persistent associations between host ancestry and mite lineages.

Microscopic mites of the genus Demodex live within the hair follicles of mammals and are ubiquitous symbionts of humans, but little molecular work has been done to understand their genetic diversity or transmission. Here we sampled mite DNA from 70 human hosts of diverse geographic ancestries and analyzed 241 sequences from the mitochondrial genome of the species Demodex folliculorum. Phylogenetic analyses recovered multiple deep lineages including a globally distributed lineage common among hosts of European ancestry and three lineages that primarily include hosts of Asian, African, and Latin American ancestry. To a great extent, the ancestral geography of hosts predicted the lineages of mites found on them; 27% of the total molecular variance segregated according to the regional ancestries of hosts. We found that D. folliculorum populations are stable on an individual over the course of years and that some Asian and African American hosts maintain specific mite lineages over the course of years or generations outside their geographic region of birth or ancestry. D. folliculorum haplotypes were much more likely to be shared within families and between spouses than between unrelated individuals, indicating that transmission requires close contact. Dating analyses indicated that D. folliculorum origins may predate modern humans. Overall, D. folliculorum evolution reflects ancient human population divergences, is consistent with an out-of-Africa dispersal hypothesis, and presents an excellent model system for further understanding the history of human movement.

Genomic signatures of geographic isolation and natural selection in coral reef fishes.

The drivers of speciation remain among the most controversial topics in evolutionary biology. Initially, Darwin emphasized natural selection as a primary mechanism of speciation, but the architects of the modern synthesis largely abandoned that view in favour of divergence by geographic isolation. The balance between selection and isolation is still at the forefront of the evolutionary debate, especially for the world's tropical oceans where biodiversity is high, but isolating barriers are few. Here, we identify the drivers of speciation in Pacific reef fishes of the genus Acanthurus by comparative genome scans of two peripheral populations that split from a large Central-West Pacific lineage at roughly the same time. Mitochondrial sequences indicate that populations in the Hawaiian Archipelago and the Marquesas Islands became isolated approximately 0.5 Ma. The Hawaiian lineage is morphologically indistinguishable from the widespread Pacific form, but the Marquesan form is recognized as a distinct species that occupies an unusual tropical ecosystem characterized by upwelling, turbidity, temperature fluctuations, algal blooms and little coral cover. An analysis of 3737 SNPs reveals a strong signal of selection at the Marquesas, with 59 loci under disruptive selection including an opsin Rh2 locus. While both the Hawaiian and Marquesan populations indicate signals of drift, the former shows a weak signal of selection that is comparable with populations in the Central-West Pacific. This contrast between closely related lineages reveals one population diverging due primarily to geographic isolation and genetic drift, and the other achieving taxonomic species status under the influence of selection.

A phylogenomic analysis of turtles.

Molecular analyses of turtle relationships have overturned prevailing morphological hypotheses and prompted the development of a new taxonomy. Here we provide the first genome-scale analysis of turtle phylogeny. We sequenced 2381 ultraconserved element (UCE) loci representing a total of 1,718,154bp of aligned sequence. Our sampling includes 32 turtle taxa representing all 14 recognized turtle families and an additional six outgroups. Maximum likelihood, Bayesian, and species tree methods produce a single resolved phylogeny. This robust phylogeny shows that proposed phylogenetic names correspond to well-supported clades, and this topology is more consistent with the temporal appearance of clades and paleobiogeography. Future studies of turtle phylogeny using fossil turtles should use this topology as a scaffold for their morphological phylogenetic analyses.

Two new species of the gorgonian inhabiting barnacle, Conopea (Crustacea, Cirripedia, Thoracica), from the Gulf of Guinea.

Two new species of Conopea (Say 1822) are described from the Gulf of Guinea: Conopea saotomensis sp. n.and Conopea fidelis sp. n. These two new species were collected from the historically isolated volcanic islands of São Tomé and Príncipe. The relationship between Conopea saotomensis sp. n., Conopea fidelis sp. n.and two other Atlantic barnacle species, Conopea calceola (Ellis 1758) and Conopea galeata (Linnaeus 1771), is examined. The methods employed are the construction of a molecular phylogeny using mitochondrial COI and nuclear H3 gene sequence data along with morphological comparisons of calcareous and cuticular body parts. It is found that Conopea saotomensis sp. n., Conopea fidelis sp. n.and Conopea calceola are most closely related to each other but the relationship among them is unresolved. Gorgonian hosts are identified. Preliminary observations show species level host specificity for Conopea fidelis sp. n.